Multiple-quantized vortices in rotating LOFF state of ultracold Fermi superfluid gas 
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A rotating ultracold S'-wave superfluid Fermi gas is considered, when the population imbalance 
5n (or equivalently the mismatch in chemical potentials Sjj.) corresponds to the Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) state in the vicinity of the Lifshitz critical point, ft is shown that under these 
conditions the critical angular velocity Qc2 in two-dimensional systems is an oscillating function of 
temperature and 5n giving rise to reentrant superfluid phases. This leads to vortex lattices with 
multiple-quantized circulation quanta. The reason for this behavior is the population by Cooper 
pairs of the Landau levels above the lowest one. 

PACS numbers: 



I. INTRODUCTION 

The Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) phase 
with a spatial variation of the superconducting order pa- 
rameter A(r) can be realized in superconductors in ex- 
ternal magnetic field (/i). The Zeeman interaction be- 
tween spins of electrons and the external field h produces 
unequal populations of the spin up and down electrons. 
This, in turn implies a mismatch in chemical potentials 
dfi = h = (/X| — ^i\)/'2, of spin-up and down compo- 
nents [Eg. 

In ultracold fermionic superfluids the pairing occurs 
between atoms in two different hyperfine states (a, 6), 
which can be loaded into a trap with unequal numbers. 
The situation is then analogous to the spin-polarized elec- 
trons in metallic superconductors. Due to the rather 
large relaxation time for transitions between the differ- 
ent hyperfine states, the numbers of both species 
and nh can be considered as fixed to a good approxi- 
mation. Since the particle losses are negligible, the pop- 
ulation imbalance Sn = Ua — rib is fixed as well. The lat- 
ter implies that the chemical potentials are mismatched; 
this mismatch is commonly characterized by the quan- 
tity 5^, = {fig^ — fJ'i,)/'2. The possibility to fix Sn means 
that the physical realization of the LOFF phase in a su- 
perfluid ultracold Fermi gas is more favorable than in 
metallic superconductors. 

In metallic S'-wave superconductors placed in an ex- 
ternal magnetic field the genuine LOFF state is usu- 
ally masked by the orbital effects due to the Lorenz 
force term in the Hamiltonian described by the mini- 
mal coupling prescription: p — > p — eA, where p and 
e are the electron momentum and charge, A is the vec- 
tor potential. To eliminate the orbital effects one needs 
specially designed materials like in a clean quasi-2D su- 
perconductor which is exposed to a magnetic field ap- 
plied parallel to the conduction planes. Another example 
arc the heterogeneous structures, where the LOFF phase 
can be realized in a much broader range of parameters 
than in bulk superconductors or superfluids. Specifically, 
in metallic superconductor-ferromagnet-superconductor 
{SFS) weak links the LOFF phase has already been re- 



alized within the ferromagnetic link. These systems are 
very promising for small-scale applications, such as the 
SQUIDs, quantum computing, etc. 

The complications mentioned above does not arise in 
superfluid two-component ultracold Fermi gases, where 
the number density is controlled in experiments, and or- 
bital effects are absent due to the neutrality of gases. 
However, the realization of the LOFF phase in super- 
fluid ultracold Fermi gases is a still unresolved issue. 
These systems are commonly placed in magnetic traps 
that create an inhomogeneous confining potential which 
plays the role of the container that keeps the ultracold 
gas in equilibrium. In such an inhomogeneous system 
a nonzero population imbalance leads to a phase sep- 
aration into a homogeneous superfluid (S) and normal 
(N) fluid. This leaves a rather narrow spatial window 
between S and N phases where the LOFF phase can 
be realized [^. The optical lattices may offer more fa- 
vorable conditions for the realization and applications of 
the LOFF phase in ultracoldgases |J]. One possible ap- 
plication, proposed in Ref. [5|, is the construction of a 
heterogeneous SNS weak link, where is an ultracold 
Fermi gas in the normal state with population imbalance 
driN ^ (and Sfij^ ^ 0), while the superfluid banks 
S are from superfluid ultracold gas (of the same atoms), 
but without population imbalance, i.e., dns = 0. It turns 
out that in such a case, depending on parameters of the 
system, one can in principle realize a 7r-weak link (a type 
of Josephson contact with the phase difference ip = n in 
the ground state), which when placed in a superfluid ring 
can produce spontaneous mass flow in the ground state 
- the so called n - SQUID [i|. Since the LOFF state 
is characterized by a negative stiffness (see below) which 
causes periodically modulated superconducting (super- 
fluid) order parameter A(r + L) = A(r), this fact may 
gives rise to a number of interesting effects. 

In this work we consider 2D systems with the LOFF 
phase in a two-component ultracold Fermi gas under ro- 
tation. The angular rotation of superfluid gases plays a 
role that is similar to the orbital effect in metallic super- 
conductors. This fact opens up the possibility to investi- 
gate in ultracold superfluid Fermi gases the interplay of 
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"orbital" and "spin" effects. In contrast to metallic su- 
perconductors these quantities can be varied in ultracold 
gases independently. Below we investigate the behav- 
ior of the critical angular velocity flc2 for the transition 
from the supcrfluid to the normal state in the vicinity 
of the Lifshitz critical point (with the temperature T* 
and the chemical potential mismatch Sfi* in the T — dfi 
plane). We show that ilc2 is an oscillatory function of 
temperature or population imbalance. These oscillations 
arise due to the occupation of (higher) Landau levels by 
Cooper pairs. The latter effect is realized if the quasiclas- 
sical condition hrtc2 < 7r(A;BTc)^//i is realized, where Tc 
is the critical temperature of superfluid phase transition 
and n is the average chemical potential. We refer to this 
effect as the quasiclassical oscillation effect. In the op- 
posite case, more precisely when hflc2 > Tc, the quantity 
ilc2 oscillates due to the population of the Landau levels 
in the normal state (the normal state Landau quantiza- 
tion) 0, i. This effect, which can be refered to as 
the quantum oscillation effect is negligible near the point 
{T*,S^j,*) and will not be studied here. 



II. ROTATING THE LOFF PHASE 

We consider a rotating Fermi gas with the angular ve- 
locity vector directed along the z-axis, fl = fiz. The 
ultracold gas is placed in a magnetic trap with the po- 
tential Vm{r) = [Afw^r^ + Mw^0^]/2, where uj and uJz 
are the trapping frequencies, M is the atomic mass. The 
BCS Hamiltonian in the rotating coordinate system, but 
expressed via coordinates of the inertial laboratory sys- 
tem, is given by 

H=^ J d'^xi:l[hoip,r,i) - ^i,]^^ 

i—a,b 



-g J d'^xi;^{x)i}f,{x)^f,{x)i}^{x), (1) 

where the summation is over the fermionic species, the 
integration involves the dimension of the space d, "^j^^ix) 
and ipai^) the fermionic creation and annihilation 
operators, g is the coupling constant, and ho is the sin- 
gle particle Hamiltonian, which will be specified below. 
The chemical potentials of the species are yit^j f, = ± Sfi, 
where ^ = (/x^ -I- fia)/"^ ^f^ = if^a ~ /^b)/2 are the 
average and the "mismatch" chemical potentials, respec- 
tively. In the following we specify the quantities fi and Sfi. 
Since the self-consistent equation for the superconduct- 
ing(superfluid) order parameter are the same for fixed 
dfj, and Sn the obtained results are applicable also to 
systems with fixed 5n. The single-particle Hamiltonian 
ho{p, r, L) in Eq. ^ is given by 

hoip,r,Lz) = ^ + V,n{r)-n-L, (2) 



where p = —ihV is the linear and L = f x p is the orbital 
momentum. If necessary /io(p,r,L) comprises also the 
optical-lattice potential Vop{r), i.e. one has 

/io(p,r,4) = ^(p-MVo)2 

+ ^MiLu' - + ^MluIz', (3) 

where Vfj = ilz x r is the velocity due to the rotation. 
In the following we consider the two-dimensional (2D) 
gas which is realized for the pancake-like trap with ujz S> 
UJ. It is also assumed that a; > $1 in order to keep the 
gas stable and quasi-homogeneous as much as possible. 
When necessary, the system can be additionally placed in 
a periodic optical-lattice potential Vop{r), which changes 
the atomic-particle spectrum from the parabolic one p^/2 
to the tight-binding like e(p), which creates favorable 
conditions for the LOFF phase Q . We study the problem 
near the second-order transition line between the normal 
(unpaired) and the LOFF phase in the T — dfi plane, 
which is reahzed for T ^ T* and Sfi ^ 5fi* [§]. It will 
be shown that near this line the critical rotation velocity 
^c2 (for the transition from the normal to the superfluid 
state) is small and the effect of rotation can be accounted 
for in the quasiclassical approximation. The latter means 
that the Landau quantization in the normal Fermi gas 
is negligible in this parameter region. At this line the 
equation for the order parameter can be linearized 

^= /dViifo(r-ri)e-*('-'-i^(P-2^^")A(r) (4) 
9 J 

with the kernel 

Xo(r-ri) =T^G,(w„,r-ri)G6(-w„,r-ri), (5) 

where G'a/f,(a;„, r — ri) are the single particle Green's 
functions, where hw^ = 77^ = 'KkBT{2m -\- 1) are 
the fermionic Matsubara frequencies. The characteris- 
tic scale of the kernel Ko{r — ri) is given by the super- 
fluid coherence length ^g. Therefore, the quasiclassical 
approximation in Eq. Q is valid if the phase change 6ip 
due to the rotation is small on the scales of the order of 
^Q. The phase change can be estimated as 

,,.2Mf'\M<n^^^ (6) 

Then, the quasiclassical condition S(p ^ 2tt gives hil ^ 
7T{kBTc)'^ / fi where /i — Ep and Ep is the (average) Fermi 
energy. When 57 is large enough so that the cyclotron 
radius i?si vp/^ of the atomic orbits is much smaller 
than ^g, i.e., when h^l > UbTc, it is necessary to take 
into account the Landau quantization of atomic motion 
in the normal state. 
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Consider next the Fourier-image of the kernel 
Xo(q) = fcsT^ / d^pGaiu;m,p)Gbi-u;m,P + ci). (7) 

In the quasi-homogeneous case, when the effect of the 
potential M(w^ — r2^)r^/2 on the single-particle spectrum 
is small, the Green's functions Ga.b are given by 



(8) 



where ^(p) — (p^/2M) — /i. After the integration over 
the energy ^ and by assuming a circular Fermi surface 
one obtains (in the following fc^ = 1) 



Xo(q,??™>0)==z7V(0) 



diLp 



2{iri^ + 5^) + qvFCOSip' 



where wj? is the Fermi velocity and A^(0) is the density of 
states at the Fermi surface. The integration in Eq. Q is 
straightforward (the contour integral runs over the con- 
tour jzj = 1 with z — exp{i</j}) and we obtain 



Ko{q) = Re 5] 



7rA^(0)T 



(10) 



Since we consider the problem near the Lifshitz critical 
point {T*,Sii*), the wave vector of the LOFF phase is 
small, i.e., q ^q^-, and it is sufhcient to make the small 
q expansion. Thus, we approximate Ko{q) « Ko{0) + 
K2q^ + K4q^, with Ko{0) = A^(0)[.9"^ - T{t,Sfl)], where 
T(t, S/2) is defined in Eq. p^ . Upon defining dimensional- 
less unit X = {q^2M)/Ec with Ec = T^/Ep, we rewrite 
this expansion as Ko{q) « Ko{0) + + k/^X'^, where 



k2 = 

^4 — 



r 



(47r)2 



Reijj 



(2) 



'27rt 



4(47r) 



rRetp 



(4) 



- + 1 
2 2nt 



and 



r(i, Sn) = Int + Re^h ( - + i — 



(11) 
(12) 

(13) 



Here, t — {T/Tco) and Sjl — Sfi/Tco, where Tco is the 
critical temperature of superfluid phase transition when 
the number densities of the two species are equal (Sn = 
0), and ipi^) = dhiT{x)/dx is di-gamma function and 
^p^^\x) = d'^ip{x)/dx^ . In case of i7 ~ one obtains from 
Eq. ([4]) the equation defining the transition line between 
the normal and the LOFF phase 

[TLitQ.Sjl) - k2E{Q) - kiE'^iQ)] = 0, (14) 

where A(r) = A, exp[Q • r] and E{Q) = 
{2Q^M)/Ec = fc2/2|fc4|. The magnitude of the LOFF 
wave vector Q is determined by the maximum of the 
function k2E{Q) + k4E^{Q), i.e., by maximizing the 



LOFF critical temperature tL{Q,6fl). The direction of 
the wave-vector is chosen spontaneously. At the Lif- 
shitz point where t* = 0.56 and Sjl* = 1.04 one obtains 
TL{t\5jj*) = and k2{t*,6fl*) = 0, where Q{t*,6fi*) = 
0. The LOFF phase is realized for t < t*and Sfl > S^i*, 
therefore near the Lifshits point one has Q(t, Sfl) ^ ^q^, 
which justifies the small-g expansion of the kernel A'o(q) 
above. 

In the case of a 2D rotating system, with fl ^ Oand 
Nil = r2 X r, there is an upper critical angular veloc- 
ity ric2(i, (^/i) below which there is a nucleation of su- 
perfluidity in the form of quantized vortices. The linear 
Ginzburg-Landau equation for the second-order normal 
state - LOFF transition reads 



Tit,Sjl)-[ -f\H 



] fr2 



A(r) = 0, (15) 



where H = {p — 2M\i-i) /2M is the Hamiltonian for 
the harmonic oscillator in the isotropic Coulomb gauge. 
To obtain flc2{t,5fl) we need the solution of the eigen- 
value problem iJA„(r) = e„A„(r). By assuming a 
rotationally infinite 2D system, such as disc, the so- 
lution of the eigenvalue problem is given by A„(r) ~ 
exp{ — (|z|/2Zo)^}, with z = x + iy and the "mag- 
netic" length Iq = y^h/4Milc2, while the eigenvalues are 
e„ — 4Slc2(?^ + 1/2), n = 0,1, 2... Note that the solutions 
A„ (r) are also eigenstates of the 2;-the component of the 
angular momentum operator = —ih{xdy — ydx) = 
h{zdz - z*dz'), i.e. 



LzA„(r) = hnAn{r). 



(16) 



The solutions with ri — 1,2... may correspond to a 
multiply-quantized vortex lattice while the case with 
n ^ corresponds to the standard Abrikosov solution 
where the number of the zeros of the order parameter and 
flux quanta in the unit cell is one. However, the calcula- 
tion of the real structure of the vortex lattice for < ric2 
requires the knowledge of the nonlinear Ginzburg-Landau 
equation and the superconducting current in the LOFF 
state, which are much more complicated than the stan- 
dard equations in the homogeneous case \^ . We expect 
that in the rotating LOFF state with n = 1,2... various 
vortex lattices may appear, where the number of zeros of 
A(r) in the unit cell can be larger than the number of 
flux quanta 11]. The solution of the eigenvalue problem 

posed by Eq. (jTS]) gives the explicit expression ri|.2''± for 
given n 



1 k2 ± yjkl - 4T|fc. 



8|/C4 



(17) 



where ujc — Ec/h. It is seen that for the fixed value 
of Sjl > Sfl* (or for fixed Sn) one has two branches 



n 



in) 
c2,+ 



> which arc functions of the temperature 

t < t* . They meet each other on the line k2{t,Sjl) — 
4:T{t, dil)\k4{t,Sjl)\ = 0, which is in fact the LOFF line 
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FIG. 1: Schematic figure of the oscillation of the critical an- 
gular velocity flc2{X; y) as the function of the variable X for 
fixed y. For the variable temperature t = T /Tco and fixed 
mismatch in chemical potentials Sjj, (or the population imbal- 
ance 5n) one has X = t and y = Sp, (or y = 5n) and vice 
versa. Xq is the LOFF critical temperature II or the critical 
mismatch Sjlj^. 

tL{tQ, Sfl) given by Eq. The lower line for the n-th 

level intersects the upper line for the n + 1-th level at the 
point ri^2^_(t„_„+i) = ^''^^\tn.n+i)- This means that 
in the temperature interval t„.„+i < t < (for fixed dfl) 
the normal state is realized for nl!^^_^^\t) < Q. < r2[.2^_(t). 
Therefore, a cascade of oscillatory (interchanging) nor- 
mal and LOFF-vortex-lattice reentrant transitions will 
arise, i.e. a sequence of transition from the normal ul- 
tracold Fermi gas to the superfluid gas with multiply- 
quantized vortex lattice. This is clearly seen in the phase 
diagram Qc2{t) shown in Fig. 1. 

We would like to stress that the phase diagram in 
Fig. 1 is also generic for the case when the temperature 
t < t* is fixed but the mismatch in chemical potentials 
Sjl > 6p,* (or Sn) is varied. In that case the variables t 
and Sjj change the roles but the oscillatory effect is the 
same. Again one has fl'"^''^{Sfl) > rtl^'^_{Sp.) and these 
two curves meet at the LOFF transition line 5jiQ{t). In 

the intervals f}^2,+^^(^A^) < < 

one has again 

the cascade of normal to the LOFF-vortex-lattice reen- 
trant transitions. 



III. DISCUSSION AND CONCLUSIONS 

Our study shows that it is possible to create a rotating 
LOFF state in superfluid ultracold Fermi gases in a quasi- 
2D magnetic trap, if the trap-frequency uj is adapted to 
be around (but still larger than) the angular velocity 



Q.c2{T,5n^) given by Eq. HT]), i.e., u > ^2^2- The lat- 
ter condition allows, first, the stability of the ultracold 
gas and second it gives rather small inhomogeneity in 
the quasiparticle spectrum. This means that the kernel 
Aro(r) in Eq. ([5]) can be calculated with the help of the 
Green's functions for the homogeneous system given by 
Eq. ([5]) . The realization of the LOFF phase is even more 
favorable if the ultracold gas is placed in a 2D optical lat- 
tice, since the existence of the van Hove singularities near 
(or at) the Fermi surface can favor the LOFF state in a 
much broader region of the phase diagram {6n,T) g. In 
the case of the rotating LOFF superfluid an additional 
weak magnetic trap, with iv > ^c2, compensates again 
the centrifugal potential due to the rotation. 

(n) 

The order of magnitude of fi^j near the Lifshitz point 

is rj^2^ - (10-1 - 10-'^)iT^/EF) and in the case of the 
superfluid Li^ where Ep ^ 2 fiK and ~ Ep/lO 

one has fl^^^ < (30 — 300) s^^. These values are in 
the range of the experimentally reached rotational fre- 
quencies realized by stirring methods for creating vor- 
tices Finally, we stress that the whole analysis was 
carried out in the framework of the BCS weak coupling 
theory. We anticipate that our analysis can be extended 
towards the Bose- Einstein condensate limit, and it should 
remain (qualitatively) valid at stronger couplings, but on 
the BCS side of the phase diagram. This expectation is 
justifled by the studies of the vortex state in metallic su- 
perconductors, for which the fermionic excitations in the 
vortex cores are realized also in the strong coupling limit 
but still on the BCS side of the phase diagram. 

In conclusion, we have studied a rotating two- 
component superfluid ultracold Fermi gas with a mis- 
match in chemical potentials of the species (or with pop- 
ulation imbalance), which is in the LOFF phase and near 
the Lifshitz point. We have shown that the critical an- 
gular velocity f2c2 in two-dimensional systems is an oscil- 
latory function of temperature or the population imbal- 
ance Sn{= Ua — rib). This effect gives rise to a cascade of 
reentrant superfluid transitions with the superfluid fea- 
turing multiple-quantized circulation quanta. The reason 
for this oscillatory effect is the population of the higher 
Landau levels {n ^ 1) by Cooper pairs in the LOFF 
state. The obtained results on the critical angular veloc- 
ity ilc2 in the LOFF state might be of interest also for 
rotating 2D metallic superconductors in the parallel mag- 
netic fleld and for color superconductors in quark mat- 
ter, where in the two-flavor color-superconducting quark 
matter the LOFF phase may compete with the gluonic 
condensate [T3| . 
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